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Completing the Calculation of BLM corrections to B — y X s j 
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Perturbative 0{a 2 s ) corrections to B(B — > X s -y) in the BLM approximation receive contributions 
from two-, three- and four-body final states. While all the two-body results are well established 
by now, the other ones have remained incomplete for several years. Here, we calculate the last 
contribution that has been missing to date, namely the one originating from interference of the 
current-current and gluonic dipole operators (K^ 130 and ). Moreover, we confirm all the 

previously known results for BLM corrections to the photon energy spectrum that involve the 
current-current operators (e.g., K^^° and i"Q? )■ Finally, we also confirm the recent findings of 
Ferroglia and Haisch on self-interference of the gluonic dipole operator (iQg °). 

PACS numbers: 12.38.Bx, 13. 20. He 



I. INTRODUCTION 

Weak radiative decay of the B meson is a well-known 
probe of physics beyond the Standard Model (SM). Cal- 
culations of its inclusive branching ratio in the SM for 
EL, > 1.6 GeV give 0,0 

B{B -> X S 7) SM = (3.15 ± 0.23) x 1CT 4 , (1) 

which agrees within 1.2a with the world average Q 

B(B -> A s7 ) cxp = (3.55 ± 0.24 ± 0.09) x 10" 4 . (2) 

The above experimental result includes a model uncer- 
tainty that is due to averaging several measurements 
with various photon energy cuts E and extrapolat- 
ing them to £q = 1.6 GeV where the theory predic- 
tion is most reliable. Measurements with energy cuts 
1.8 GeV < E a < 2.0 GeV have significantly 

smaller background-subtraction errors than those with 
E Q = 1.7 GeV [|. More work at balancing model- 
dependence and background-subtraction uncertainties is 
necessary in the future to obtain accurate experimental 
averages. 

As far as the SM calculations are concerned, further 
improvements require another critical rc-analysis of non- 
perturbative effects 0, as well as a full perturbative 
0(a 2 s ) evaluation of T(b — > Xfj), where Xf stands for 
s, sg and sqq partonic states (q = u,d,s). Such cal- 
culations are most conveniently performed in the frame- 
work of an effective low-energy theory that arises from 
the SM via decoupling of the W boson and all the 
heavier particles. So long as higher-order electroweak 
and/or CKM-supprcsscd effects are neglected, the rele- 
vant flavor-changing weak interactions at the renormal- 
ization scale fib ~ wif,/2 are given by 



4G 



-'weak 



■a 



i=l 



Ci(fj,b)Qi, 



(3) 



where Qi denote either dipole-type or four-quark opera- 
tors (see below), and Ci(fib) stand for their Wilson coef- 
ficients. 



Following Refs. [1, we shall normalize the radia- 
tive decay rate to the charmless semilcptonic one, and 
parametrize their rato in terms of symmetric matrices 
Kij(fib, Eq) as follows: 



T[b -y X u ev] 



v t *v tb 



2 a 8 



i 1 J n ' • • 



(4) 



where Cf s are certain linear combinations of Cj, see 
Eq. (5) of Rcf. 0. Evaluation of all the Cf{y, b ) up 
to the Ncxt-to-Next-to-Lcading Order (NNLOJin QCD 
has been already completed several years ago [10] ■ 
In the perturbative expansion of Kij 



Ki 



K 



(o) 



4tt » 



(2) 



(5) 



all the 0(1) and G(a s ) terms are known since a long 
time [Til] - As far as are concerned, the so-called pen- 
guin four-quark operators Q3, . . . , Qe can be neglected 
thanks to smallness of their Wilson coefficients. We can 
restrict our attention to 



(s L7fl T a c L )(c L ^T a b L ), 
(sL7M c i)(ci7'' & i), 

P 

mbisLcr^b^F^, 
m b {s L a^T a b R )G%, 



Q7 

i.e. consider K, 



16tt 2 
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(2) 



16tt 2 

with i,j e {1,2,7 



(6) 



only. 



At present, are known in a complete manner 112- 
[l| for (ij] = (77) and (78), while the other cases are 
estimated [15|-[l3] using the BLM [lj| approximation. In 
Ref. (H, non-BLM contributions to the decay rate have 
been calculated in the m c 3> rribjl limit, and then inter- 
polated downwards in m c assuming that they vanish at 
m c = 0. Such a treatment of non-BLM NNLO correc- 
tions in the evaluation of Eq. ([1} still remains the cur- 
rent state-of-art for the numerically important but yet 
unknown and . 
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FIG. 1: Feynman diagrams that determine A' 



(2)/3 



TABLE I: Present status of K, 



(2)/3 



calculations. 



II. CALCULATION OF K 



AND A 



(2)/3 



(2) 

The BLM and non-BLM contributions to K) 4 are de- 

noted by K^^ and AT^ 2 ' rGm , respectively. The latter are 
independent on ni (the number of massless quark flavors) , 
while the former are proportional to /?o = 11 — |(n; +2). 
In practice n; = 3 because masses of the light q = u,d, s 
quarks are neglected in loops on the gluon lines in b — > sj 
and b — > sgj, as well as for external qq pairs in b — > 
s 9*l S QQ1- Although masses of the c and b quarks 
are not neglected, all the quantities in Eq. (|5|) are MS- 
renormalized at in the five-flavor theory, which justi- 
fies the use of five-flavor /3 in A^ 2)A) . Effects of non-zero 
values of m c and mb in loops on the gluon lines are known 
from Rcfs. [Hd] for all the K\f with i, j g {1, 2, 7, 8}. 
No real cc pair production is included in b — > X^j by 
definition, while bb production is kinematically forbidden 
anyway. 

Contributions to from the b — > s"/ channel arise 

for (ij) = (17), (27), (77) and (78) only. They were orig- 
inally calculated in Ref. [l6[ ■ Three- and four-body final 
state contributions (6 — > sgj and b — > sg*j — > sqqj) for 

all the i,j € {1, 2, 7} cases and for K^^° were evaluated 



first in Ref. [L_ 
Ferroglia and Haisch [ttI ]. 



Recently, K^f 00 



has been found by 



Determination of K^ 00 and K^^° amounts to eval- 



,(2) 



and K$. The latter 



uating n;-dependent parts of K ls 
originates from interference of decay amplitudes gener 
ated by the current-current operator Qi and the gluonic 
dipole operator Qs- The contributing Feynman diagrams 
are most conveniently presented using Cutkosky rules (2p| 
as four-loop propagator diagrams with unitarity cuts. 
They are displayed in Fig. [1] In dimensional regulariza- 
tion, no diagrams with cuts through the gluon lines need 
to be considered because the massless gg-loop integral is 
scaleless for an on-shell gluon, which implies that all such 
diagrams vanish. If Qi is replaced by Q±, the color fac 
tor gets modified according to T a T b T a T b 
which leads to a simple relation 



6 J ' 



A' 



(2)/?o 
18 



-A". 



(2)/3o 
28 



(7) 



Following the conventions introduced in Ref. [15(, we 
exclude the diagrams depicted in Fig.[2]from the BLM ap- 
proximation despite their n;-dcpcndcncc. They are corre- 
lated via renormalization group with tree-level b — > sqqj 
matrix elements of the neglected four-quark operators 
Q3, Qq. Excluding those diagrams from the BLM cal- 
culation is indeed reasonable. No other n;-dependent dia- 
grams arise because the (32-generatcd charm loops vanish 
if the on shell photon alone is emitted from them. 



In the present paper, we provide the last two missing 



contributions, namely AT^g^ 



and K, 



(2)A) 
28 



Moreover, we 



(11), (12), (22), (17) and 



confirm the results for (ij) 

(27) from Ref. [F 

Table |T] summarizes the present status of K 1 ^ 00 
tions. 



as well as for K^ 00 from Ref. [Li 



calcula- 



The article is organized as follows. In Sec. [TTI our 
evaluation of A^g^ and AT^g^ is presented. Sec. IIIII is 
devoted to the remaining contributions that involve the 
current-current operators (Q\ and Q2)- Self-interference 
of the gluonic dipole operator Qg is considered in Sec. lIVI 
We conclude in Sec. IVl 






FIG. 2: Feynman diagrams excluded from K. 



(2)ft> 
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In our actual evaluation of -K28 > the method of Smith 
and Voloshin [2l| has been applied. It amounts to consid- 
ering lower-order diagrams that are obtained from those 
in Fig. [T] by removing the qq loop from the gluon propa- 
gators. However, an arbitrary auxiliary mass of the gluon 
needs to be introduced. Next, integration over the gluon 
mass should be performed according to the formulae of 
Ref. [2l|. 

We have carried out the calculation via direct integra- 
tion over the 3-body partly massive phase space (m s = 0) 
that is conveniently parametrized in terms of two vari- 
ables: u = 2{p b p 1 )/m 2 and s = (p g + p~ f ) 2 /m 2 . Explicit 
results from Sec. 4 of Ref. [22| for the one-loop Qi ampli- 
tude with an external off-shell gluon have appeared to be 
useful. Once the Dirac algebra is performed, we are left 
with precisely the same two Feynman parameter integrals 
as in Eqs. (4.21) and (4.22) of that paper, namely 



F b (s,z,v) 



F g (s,z,v) 



1 pi 



dx / dy — — 

Jo v + ys- 

1 - 1 -xy 



dx I dy — 

Jo v + ys 



(8) 



o(l-x) 



TO giuon l m l ■ Considering 



where z = m 2 /m 2 and v 

them here in D = 4 dimensions is sufficient because the 
calculation is free of ultraviolet, infrared or collinear di- 
vergences. Integrations over the two Feynman param- 
eters x and y and over the phase-space variable u are 
performed analytically in a straightforward manner. The 
remaining two integrations (over s and v) have been com- 
pleted numerically. More details will be presented else- 
where [23j . 

For K\f Po with i,j ^ 7 we shall use the following 
notation (consistent with Ref. 0): 



K 



(2)A> 



= 2(1 + <%)A) <$(S)L h + h%>(S) , (9) 



(2) 



where 8 = 1- 2E Q /m b , L b = In (// 2 /mf), and 4>\j\S) 
are the well-known NLO brcmsstrahlung functions col- 
lected in Appendix E of Ref. [§[ . 

Our final result for the function h y 2S ' (S) reads 

h^s(S) = 0.02605 + 0.1679 (5- 0.1970 5 2 

+ (-0.03801 + 0.6017 <5- 0.7558 S 2 ) z^ 



(2.755- 10.03 5+ 11. 27<5 2 ), 
(-27.05 + 68.47(5- 72.51 5 2 



z 2 



+ (85.87 - 289.3 5+ 297.7 5 2 ) z 2 
+ (-91.53 + 399.8 S- 399.9 5 2 ) z? . 



(10) 



The above expression is a numerical fit that remains ac- 
curate in the ranges < z < 0.13 and 0.2 < S < 0.6. 
These ranges will also be valid for the fits in Sec IIHI The 
central values used in Eq. (TTJ) arc 5=1 — 2(1.6/4.68) ~ 
0.316 and z = [m c (1.5GcV)/4.68] 2 ~ 0.0584. 



Eq. (|10p is the main new result of the present pa- 
per. Its numerical effect on the branching ratio turns out 
to be minisculc (below 0.1%). However, the purpose of 
the present calculation is not finding sizeable effects but 
rather removing several minor uncertainties that had to 
be taken into account in Refs. 0, Q in estimating the 
±3% perturbative error that was unrelated to the m c - 
interpolation. 



III. OTHER CONTRIBUTIONS FROM 
CURRENT-CURRENT OPERATORS 

Let us now consider k\ 2)0 ° for (ij) e {(11), (12), (22)}. 
The three Feynman diagrams to be calculated in the (22) 
case are obtained from the left parts of the cut diagrams 
in Fig. [1] by forming all the possible interference terms. 
The cases (11) and (12) differ from (22) by color factors 
(analogously to Eq. ([7])), namely 



K. 



(2)/3o 



-6i^2 )Al =36i^ 2)A) . 



(11) 



As before, only the diagrams with both the photon and 
the gluon coupled to the charm loop are included in the 
BLM approximation for the b — > s<?*7 — > sqq^f channel. 

Using precisely the same methods as in Sec. UH we 
obtain the following numerical fit: 



0.01370 + 0.3357 5 - 0.08668 S 2 
+ (0.3575 + 1.825(5- 0.3743 6 2 ) z? 
+ (-2.306- 5.800 5- 6.226 S 2 ) z 
+ (3.449 -0.5480 <5+ 17.27 <5 2 ) zi (12) 

Similarly, for the photonic dipole (Q7) and the current- 
current operator interferences, we find 

hfh 5 ) = -0.1755 - 1.455(5+ 1.119 (5 2 

+ (0.7260- 7.230 (5 + 5.977 <5 2 )z5 
+ (13.79 + 113.7(5- 100.4 <5 2 )z 

+ (-145.1 - 307.1 <5 + 388.5 5 2 ) z% 

+ (475.2 + 313.0(5- 775.8 <5 2 )z 2 

+ (-509.7- 126.1 <5 + 646.2 (5 2 )zi, (13) 

together with K^ 00 = -\K {2 ^° . However, the two- 
body contribution T is non-vanishing in this case, so in- 
stead of Eq. ([9]) one has 



K$ Po = T 



2c/>. 



(2)/3o _ 



= T + 2/3 



7 (2) 
'27 



(14) 



Explicit formulae for T can be found in Ref. Q ■ 

To compare our expressions in Eqs. (TT2"|) and (TIB")) to 
Ref. [lU , we made use of their results provided to us in 
the form of a numerical grid [24| . The grid described 
contributions to the differential photon energy spectrum 
in the ranges 0.2 m b < E 7 < m b /2 and < z < 0.13. 



4 



The originally published fits in Eq. (12) of Rcf. [l5j were 
valid in more narrow ranges, especially in the case of z. 
After performing an accurate fit to the grid, we integrated 
the spectrum over from Eq to TOb/2, and then com- 

(2) (2) 

pared the outcome to our results for h 22 '(8) and h 27 '(S) 
in the ranges 0.2 < S < 0.6 and < z < 0.13. A perfect 
agreement was found immediately at the first attempt to 
perform such a comparison, indica ting undoubtedly that 
we can confirm the results of Ref . [15| . 

(2) 

It is interesting to observe that h 22 affects the branch- 
ing ratio by +1.9%, which remains within the assumed 
±3% uncertainty for all such effects in Eq. (J]). Neither 

(2) (2) 

this contribution nor the smaller ones from h 27 and h 7S 
were included in the analysis of Refs. [l], H|. 



IV. GLUONIC DIPOLE OPERATOR 
SELF-INTERFERENCE 



spectrum is found by differentiating K S8 with respect 
to 6. Doing so, we find perfect agreement with the 
very recent article of Ferroglia and Haisch [ItJ ■ An ex- 
tended discussion of collinear divergences can be found 
there, which adds new elements to the previous analy- 
ses in Refs. [3, Replacing the perturbative collinear 
regulator m s by a physical hadronic one can hardly 
be performed in a quantitatively precise manner given 
our poor knowledge of the QCD bound state properties. 
Fortunately, the gluonic dipolc operator self-interference 
undergoes significant suppression in B — > X s j due to 
{QdCs/Cj) 2 ~ 1/36, as well as the relatively high photon 
energy cut Eq ~ m;,/3. Evaluation of K^ 13 " provides 
just a check that no unexpected large numerical factors 
overcome this suppression at the perturbative level. The 
overall effect of K^^° on the b —> Xfj decay width does 
not exceed 0.2%. 



The three Feynman diagrams that matter for K^^° 
are obtained from the right parts of the cut diagrams 
in Fig. [T] by forming all the possible interference terms. 
Contrary to the previously discussed cases, a collinear 
divergence arises here, and a non- vanishing mass of the s- 
quark must be retained in the line from which the photon 

is emitted. In analogy to the NLO calculation of -Kgg , we 
shall keep this mass whenever it can produce ln(mb/m s ), 
but neglect all the power corrections (m s /mb) n . 

Note that only b —> sg*"f sqq^ matters for K 1 ^^ 1 
which means that no photon emission from the qq pair 
needs to be considered. Amplitudes with such an emis- 
sion would not be proportional to ni but rather weighted 
with the quark electric charges. Consequently, the quarks 
into which the gluon fragments can be kept massless from 
the outset, even if they are the s quarks. 

As the interfering amplitudes are tree-level here, all 
the phase-space integrals and the ficticious gluon mass 
integral can be performed analytically. We obtain 

h$(6) = ^{[(l + i£)£m<5-6m(l-£)-2Li 2 (l-5) 

' ' f'S 

+ (5-2 In J) [Lia(l-<f)- - ^tt 2 6 (2+5) 
+ [±6 + \5 2 - Hl-5)} In 2 5+(l§- i^ 2 ) x 
x ln(l- ( 5) + (-f|-i|<5+|(5 2 ) ( 51n(5 
+ 3% S +%P-%P}. (15) 
The corresponding contribution to the photon energy 



V. CONCLUSIONS 

The NNLO QCD corrections to 6 -4 Xfj in the 
BLM approximation receive contributions from b — > S7, 
b — > sg7 and b — > s.g*7 — > sqqj. The former results are 
well established by now, while the b —> sgj BLM ampli- 
tudes vanish in dimensional regularization. In this arti- 
cle, we have calculated the last missing b — > sg*j — > sqqj 
contributions, namely K^^° and K^^ . In addition, 
we have confirmed all the previously known results for 
BLM corrections to the photon energy spectrum that in- 
volve the current-current operators, as well as the re- 
cently found ifgg^ . Numerical effects of all these quan- 
tities on the branching ratio remain within the ±3% per- 
turbative uncertainty estimated in Refs. [HQ- 
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